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1. Introduction 

One of the hardest problems in lattice QCD has been the simulation of QCD with a nonzero 
baryon chemical potential pL. Some progress has been made recently in lattice simulations with 
small chemical potentials at temperatures around the chiral symmetry restoration transition [fit], 
however for large /I the main progress has come from analytic approaches. Much of the analytic 
work at large /J. has focused on the phase diagram for /J. larger than the critical pL c for chiral sym- 
metry restoration [g]. Here we will instead focus on the low-energy properties of QCD for fx < [i c 
where chiral effective theories can provide exact results. 

The low-energy eigenvalue density of quenched QCD with a chemical potential was first de- 
rived from chiral effective partition functions using the replica limit of the Toda lattice equation 
[g]. The same density was later obtained by exact diagonalization of a chiral random matrix theory 
(RMT) using the orthogonal polynomial method The agreement between the results that were 
derived from two different models suggests that there is a common universality class to which both 
belong. Comparisons to the quenched eigenvalue density for lattice QCD with a chemical potential 
were carried out previously [Q] using a slightly different theoretical result [^] that turned out to 
be not quite the correct result for the chiral universality class. A later comparison of the correct 
density to lattice QCD data found very good agreement [0]. 

The solution of the RMT also provided results for all eigenvalue correlations and was extended 
to the unquenched case. The differences between the quenched and unquenched results have al- 
ready provided much insight into the nature of chiral symmetry breaking in finite density QCD [j|]. 
Here we look only at quenched eigenvalues since they are easily obtained and save the comparison 
of unquenched eigenvalues for a future work. We note that recently comparisons of unquenched 
Dirac eigenvalues for two-color QCD with lattice results were made For three-color QCD, 
all the low-energy eigenvalue correlations in a finite volume V with topological charge v can be 
obtained from the kernel 

lxvl v+1 / 7 \x\ 2 \ / ivl 2 \ Mi+h r 1 

= A$wrr\wr\wr^ L <>■» 

with x = xLV, y = yLV and jx = jxF^/V [|J]. The low-energy constants £ and F are the tree level 
chiral condensate and pion decay constant that appear in the chiral Lagrangian. The density is then 
p(z) = J^(z,z) and the two-point correlation function used below is Y2(x,y) = \ J^(x,y)\ 2 . 

Since ( [ET| ) is obtained from the zero-momentum limit of the chiral Lagrangian, it is only 
valid as long as the higher momentum modes can be neglected. At pL = this approximation 
is valid as long as the eigenvalues are smaller than the scale E c ~ F 2 fL\^V [[RJ which is the 
equivalent of the Thouless energy in condensed matter systems. The average eigenvalue spacing 
at the origin for jj. = is known to be A = % /LV [JTTJ] so that the dimensionless Thouless energy 
(conductance) is g = E c /A ~ F 2 \/V (ignoring numerical factors). Here we examine how well 
the analytic predictions for the eigenvalue correlations agree with quenched lattice simulations, 
identify the Thouless energy and see how it varies with jU and V. 

2. Determining the low-energy constants by fitting the density 

The lattices used for this study were all generated using the standard plaquette action at lattice 
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Figure 1: Various visualizations of the eigenvalue density for quenched QCD on a 6 4 lattice for different /i. 
The fits are from the analytic prediction (1.1). The fits for other lattice sizes and jj, were also very good. 



coupling j8 =5.0. The eigenvalues were obtained from the staggered Dirac operator with different 
chemical potentials. The lattice sizes and values of pL used along with the number of configurations 
studied are given in table ||. The smallest value of }X for each volume was chosen to keep the value 
of jJ-VV that appears in the analytic result (LI) constant. More details on the choice of parameters 
and calculation of the eigenvalues can be found in []|]. 

First we need to determine the parameters £ and F that appear in the spectral correlations. 
This is easily done by fitting to the eigenvalue density. One approach is to fit the density along thin 
strips in the complex plane that are essentially one dimensional as was done in []|]. In order to use 
as much data as possible we take a different approach. We fit to the integrated density 



l(r,0) 



s ds d(j> p(z : 
o Jo 



se 



(2.1) 



size 


6 4 


6 4 


6 4 


6 4 


6 4 


6 4 


6 4 


6 4 


8 4 


8 4 


to 4 


10 4 


jla 


0.006 


0.008 


0.03 


0.05 


0.1 


0.2 


0.4 


1.0 


0.003375 


0.2 


0.00216 


0.2 


configs 


15500 


508 


19000 


15802 


15346 


18500 


15346 


15342 


20000 


18609 


4000 


4000 


&'max 


0.004 


0.004 


0.004 


0.004 


0.008 


0.015 






0.001 


0.005 


0.0004 


0.001 



Table 1: Lattice simulation and fitting parameters. 
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Figure 2: Fitted values of the low-energy constants E and F for the different values of jJ, and V. All the 
fitted values are consistent with statistical errors between 0.3-3%. 



To avoid complications involved with fitting to two variables we look at two different one dimen- 
sional functions. First is the radial dependence /(r, n) which reduces to just the integrated density 
along the imaginary axis when pL = 0. Here we only need to integrate to % due to chiral symmetry. 
We also look at the angular dependence given by I{r max , 6) with r max fixed. We will compare the 
lattice data to the analytic expression at v = since the staggered fermions at large lattice spacing 
do not see other topological sectors. 

In figure [I] we show a few examples of the fits for various }JL. In all cases the fits were obtained 
using the integrated density just described, though for \ia = 0.006 we plot the density given by 
either dl(r max , 6)/dd or dl(r, n)/dr for better illustration. The values of r max we used in the angular 
dependence (see table 0) are also the maximum values of r we used in the radial fits. This was 
determined by increasing r max until the % 2 /DOF of the radial fit was around one. 

At /» = the eigenvalues are purely imaginary and the angular dependence becomes a delta 
function at 6 = n/2. The radial dependence clearly shows the characteristic oscillations of the 
chiral RMT ensemble. For small pL the angular dependence moves from a delta function to a 
sharply peaked Gaussian while there is little change in the radial dependence. As /I is increased 
more, the angular dependence becomes broader and approaches a constant. Meanwhile the radial 
dependence loses the oscillations as }X is increased. 

Since we are using staggered fermions we actually have four "tastes" in the simulation instead 



of just one flavor. To correct for this we can simply replace V in ( |L1[ ) with 4V []12[]. In that paper 
the authors also demonstrate how to obtain values for £ and F, however they use the eigenvalue 
correlation function for imaginary isospin chemical potential. In figure || we show the obtained fit 
values for the different volumes and values of }l used. We cannot obtain independent fits for £ and 
F for all values of jJ. because in the large-/t limit the density reduces to a function only of the ratio 
L/F All the values of £ and F that we were able to obtain are consistent with each other. 



3. Finding the Thouless energy by examining the eigenvalue correlations 



Now that the parameters have been fixed, we look at the eigenvalue correlations to see what 
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Figure 3: Number variance E2 (f) versus the average number of eigenvalues r in a circle of varying size 
around the origin. The dashed lines are the analytic curves using the values of the low-energy constants 
obtained from the density fits. 

the effective range of the zero-momentum theory is. One convenient quantity to use is the number 
variance. For complex eigenvalues this can be denned for some region A of the complex plane as 



E 2 (r) = r- / d 2 zi [ d 2 z 2 Y 2 (zi,Z2) 

J A J A 



(3.1) 



with r = J A d 2 zp(z) fll3|]. Once the values of £ and F are set there are no more free parameters. 
Here we take A to be circles centered on the origin with varying radius. 

In figure || we show the number variance versus level number (r above) on 6 4 lattices over a 
range of values of /J.. For all values of \la < 0.4 we see very good agreement between the lattice 
results and the theoretical prediction up to the scale of a few eigenvalues. For smaller jU distinctive 
oscillations are clearly seen in both the lattice data and theoretical result. As /I increases the 
oscillations go away and the number variance approaches a curve given by taking jx to infinity in 
(LI). This can be seen in the lower right plot of figure || for values of \ia = 0.1, 0.2 and 0.4 which 
all agree for small values of r and eventually turn away. However for \ia = 1 the number variance 
jumps up dramatically and does not agree with the theory. This is most likely due to saturation of 
the particle number on the finite lattice. Here we expect the eigenvalue statistics to move from RMT 



to Poisson as was seen in [|14|]. By estimating the Thouless energy as the point where the number 
variance on the lattice begins to deviate from the theoretical curve, we can see that it increases with 
jU from about 3 to about 10 eigenvalues. The precise form of the \i -dependence needs further study. 

Next we look at the scaling of the Thouless energy with V. In figure ^ we show the number 
variance for the smallest values of pL for each volume. Here we clearly see that the Thouless energy 
is increasing with the volume. We do not have enough statistics at the largest volume to accurately 
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Figure 4: Number variance for small /xaanddif- Figure 5: Number variance for large jxa = 0.2 
ferent V. and different V. 

measure the exact scaling, but it does seem to be consistent with the \/V behavior expected at \i = 0. 
For larger \la = 0.2 (see figure ||), the Thouless energy also appears to increase with volume but 
again we do not yet have enough statistics for an accurate determination of the scaling. 

4. Conclusions 

We have shown that the low-energy eigenvalue density and two-point correlations of quenched 
QCD with a chemical potential are accurately determined from the chiral effective theory. From 
fitting the density to lattice calculations we easily obtained the low-energy constants £ and F. We 
also could clearly see the Thouless energy and its general dependence on V and jj,. In future work, 
this dependence should be quantified. We should also move to weaker coupling and larger lattices. 

We would like to thank JJ.M. Verbaarschot for helpful discussions that inspired this work. 
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